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Abstract — This paper considers the problem of robust stabil- 
ity for a class of uncertain nonlinear quantum systems subject 
to unknown perturbations in the system Hamiltonian. The 
case of a nominal linear quantum system is considered with 
non-quadratic perturbations to the system Hamiltonian. The 
paper extends recent results on the robust stability of nonlinear 
quantum systems to allow for non-quadratic perturbations 
to the Hamiltonian which depend on multiple parameters. A 
robust stability condition is given in terms of a strict bounded 
real condition. This result is then applied to the robust stability 
analysis of a nonlinear quantum system which is a model of an 
optical parametric amplifier. 

I. Introduction 

In recent years, a number of papers have considered the 
feedback control of systems whose dynamics are governed 
by the laws of quantum mechanics rather than classical 
mechanics; e.g., see [1] — [13]. In particular, the papers [10], 
[14] consider a framework of quantum systems defined in 
terms of a triple (S,L,H) where S is a scattering matrix, 
L is a vector of coupling operators and H is a Hamiltonian 
operator. 

The papers [15], [16] consider the problem of absolute 
stability for a quantum system defined in terms of a triple 
(S,L,H) in which the quantum system Hamiltonian is 
decomposed as H = Hi + H2 where Hi is a known nominal 
Hamiltonian and H2 is a perturbation Hamiltonian, which is 
contained in a specified set of Hamiltonians W. In particular 
the papers [15], [16] consider the case in which the nominal 
Hamiltonian Hi is a quadratic function of annihilation and 
creation operators and the coupling operator vector is a linear 
function of annihilation and creation operators. This case 
corresponds to a nominal linear quantum system; e.g., see 
[4], [5], [7], [8], [13]. The results in [15], [16] have recently 
been extended to allow for uncertainties in the coupling 
operator L [17]. Also, the results of [15] have been applied 
to the robust stability analysis of a quantum system which 
consists of a Josephson junction in a resonant cavity [18]. 
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In the paper [15], it is assumed that H2 is contained in a set 
of non-quadratic perturbation Hamiltonians corresponding to 
a sector bound on the nonlinearity. In this case, [15] obtains 
a robust stability result in terms of a frequency domain 
condition. Also, the paper [16] restricts attention to quadratic 
perturbation Hamiltonians. In this case, which corresponds 
to linear perturbed quantum systems, a frequency domain 
robust stability condition is also obtained. 

An example considered in the paper [16] involves the 
robust stability analysis of a quantum system consisting of a 
linearized optical parametric amplifier (OPA). Optical para- 
metric amplifiers are widely used in the field of experimental 
quantum optics; e.g., see [19]— [21]. In particular, they can 
be used as optical squeezers which produce squeezed light 
which has a smaller noise variance in one quadrature than 
the standard quantum limit. This is at the expense of a 
larger noise variance in the other quadrature; e.g., see [19]- 
[24]. Such an OPA can be produced by enclosing a second- 
order nonlinear optical medium in an optical cavity; e.g., 
see [20], [22]-[24]. Thus, an OPA is an inherently nonlinear 
quantum system. However, the paper [16] only dealt with 
linear perturbed quantum systems and thus the results of this 
paper could only be used to analyze the robust stability of a 
linearized version of the OPA. Furthermore the results of [15] 
on nonlinear perturbed quantum systems cannot be directly 
applied to the OPA system since the results of [15] only deal 
with scalar nonlinearities but the nonlinearity in the OPA 
model is dependent on two variables; e.g., see [20], [21]. In 
this paper, we extend the result of [15] on the robust stability 
of nonlinear quantum systems to allow for non-quadratic 
perturbations in the Hamiltonian which depend on multiple 
variables. This enables us to analyze the robust stability of 
the OPA nonlinear quantum system. 

II. Robust Stability of Uncertain Nonlinear 
Quantum Systems 

In this section, we describe the general class of quantum 
systems under consideration. As in the papers [10], [14]- 
[16], [25], we consider uncertain nonlinear open quantum 
systems defined by parameters (S, L, H) where S is the 
scattering matrix which is typically chosen as the identity 
matrix, L is the coupling operator and H is the system 
Hamiltonian operator which is assumed to be of the form 
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Here a is a vector of annihilation operators on the underlying 
Hilbert space and a# is the corresponding vector of creation 



operators. Also, M G 
form 
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-<2rix2n 



is a Hermitian matrix of the 
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and Mj = M[, M 2 = M 2 . In the case vectors of operators, 
the notation t refers to the transpose of the vector of adjoint 
operators and in the case of matrices, this notation refers 
to the complex conjugate transpose of a matrix. In the case 
vectors of operators, the notation * refers to the vector of 
adjoint operators and in the case of complex matrices, this 
notation refers to the complex conjugate matrix. Also, the 
notation * denotes the adjoint of an operator. The matrix M 
is assumed to be known and defines the nominal quadratic 
part of the system Hamiltonian. Furthermore, we assume 
the uncertain non-quadratic part of the system Hamiltonian 
/(z, z*) is defined by a formal power series of the form 
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[ z\ Z2 ... z m ] is a vector of operators on the un- 
derlying Hilbert space defined by 
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Also, we write 
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It follows from this definition of f(z, z*) that 
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and thus f(z, z*) is a self-adjoint operator. The term f(z, z*) 
is referred to as the perturbation Hamiltonian. It is assumed 
to be unknown but is contained within a known set which 
will be defined below. 

We assume the coupling operator L is known and is of 
the form 
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The annihilation and creation operators are assumed to 
satisfy the canonical commutation relations: 
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where J 
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■ e.g., see [6], [11], [13]. 

To define the set of allowable perturbation Hamiltonians 
/(•), we first define the following formal partial derivatives: 
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and for given constants 7 > 0, Si > 0, S 2 > 0, we consider 
the sector bound condition 



E 



df(z,z*)*df(z,z*) 



dzi 



i=l 

and the condition 

p 



dz 



1 p 

< — E z ^ + Si 



y d 2 f( Z ,Z*) * d 2 f(z,Z*) 



(9) 



(10) 



i=l 



Then we define the set of perturbation Hamiltonians VV as 
follows: 
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/(•) of the form (0 such that 
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conditions ((9) and ( TTOb are satisfied 

Note that the condition ( TTOb effectively amounts to a global 
Lipschitz condition on the quantum nonlinearity. 

As in [15], [16], we will consider the following notion of 
robust mean square stability. 

Definition 1: An uncertain open quantum system defined 
by (S, L, H) where H is of the form (Q), /(•) € W, and L 
of the form (O is said to be robustly mean square stable if 
there exist constants c\ > 0, c 2 > and C3 > such that 
for any /(■) G W, 
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denotes the Heisenberg evolution of the 



vector of operators # ; e.g., see [14]. 

We will show that the following small gain condition is 
sufficient for the robust mean square stability of the nonlinear 
quantum system under consideration when /(■) G W: 



1) The matrix 

F = 

2) 



-iJM JN^JN is Hurwitz; 

2 



E*Y< (si ~ Fy 1 JZE 1 



< 



(13) 



(14) 



where £ = 



/ 
/ 

This leads to the following theorem. 

Theorem 1: Consider an uncertain open nonlinear quan- 
tum system defined by (S, L, H) such that H is of the form 
©, L is of the form (O and /(•) € W. Furthermore, assume 
that the strict bounded real condition (I131 l. (TBI is satisfied. 
Then the uncertain quantum system is robustly mean square 
stable. 

In order to prove this theorem, we require the following 
definitions and lemmas. 

Lemma 1 (See Lemma 3.4 of [14].): Consider an open 
quantum system defined by (S, L, H) and suppose there ex- 
ists a non-negative self-adjoint operator V on the underlying 
Hilbert space such that 

-i[V,H] + ^[V,L] + ^[L\V}L + cV < A (15) 

where c > and A are real numbers. Then for any plant 
state, we have 

(V{t)) < e~ ct (V) + -, V< > 0. 

c 

In the above lemma, [-, •] denotes the commutator between 
two operators. In the case of a commutator between a scalar 
operator and a vector of operators, this notation denotes the 
corresponding vector of commutator operators. Also, V(t) 
denotes the Heisenberg evolution of the operator V and (■) 
denotes quantum expectation; e.g., see [14]. 

We will consider quadratic "Lyapunov" operators V of the 
form 
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where P G c2nx2n is a positive-definite Hermitian matrix 
of the form 
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Hence, we consider a set of non-negative self-adjoint oper- 
ators V defined as 



V 



V of the form ( fT6l l such that P > is a 
Hermitian matrix of the form ([PTt 
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Lemma 2: Given any V G V, then 

IM = [zi, [zi, V}} = [z*, [z*,V]}* = -EiSJPJEf, (19) 

which are constants for i = 1, 2, . . . ,p. 

Proof: The proof of this result follows via a straightforward 

but tedious calculation using ((6]). □ 



Lemma 3: Given any V G V, then 
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and the constants m are defined as in ( fT9l . 

Proof: First, we note that given any V E V, i E {1,2, . . . ,p}, 

and k > 1, 
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Also using Lemma |2] it follows that for any n > 1, 
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Therefore using (l22l and (l23l . it follows that 
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which holds for any i G {1, 2, . . . ,p] and fc > 0. Similarly 
for any j G {1, 2, . . . ,p} and I > 0, 
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Now it follows from (fJTJ that condition d20b is satisfied, 
respectively. □ 
Lemma 4: Given V € V and L defined as in (|5j, then 
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Proof: The proof of these identities follows via straightfor- 
ward but tedious calculations using ([6}. □ 
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Hence using Lemma [4] we obtain 
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where F = -iJM - \JNUN. 

We now observe that using the strict bounded real lemma, 
([PTb and (fT4l imply that the matrix inequality 

F^P + PF + APJY,E T E*Y 1 JP+^rY>E T E*Y l < 0. (29) 

will have a solution P > of the form illh e.g., see [8], 
[26]. This matrix P defines a corresponding operator V E V\ 
as in (TT6b . From this, it follows using (f28b that there exists 
a constant c > such that 
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Substituting ([32]), ||33}, (O and (O into ([31}, it follows that 
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Furthermore, [V, Zj]* = z*V — Vz* = [z* ,V] since V is 
self-adjoint. Therefore, 
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Note that the strict bounded real condition (fl4] i can be 
simplified according the following lemma. 

Lemma 5: The strict bounded real condition (fT4t is sat- 
isfied if and only if the following equivalent strict bounded 
real condition is satisfied: 
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Proof: First note that since £ = £ , we can write 
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Therefore, it follows from (139) that 
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From this it immediately follows that the condition ( |38| > is 
equivalent to the condition < TT~4T >. This completes the proof of 
the lemma. □ 

III. Robust Stability Analysis of an Optical 
Parametric Amplifier System 

In this section, the nonlinear quantum system under con- 
sideration is the model of an OPA. An OPA consists of 
an second-order nonlinear optical medium enclosed in an 
optical cavity. The second order nonlinear optical medium is 
referred to as a x^ medium and allows for coupling between 
a fundamental electromagnetic field and a second harmonic 
electromagnetic field; e.g., see [20], [21]. The construction 
of an OPA is illustrated in Figure Q] A standard (S,L,H) 



Partially reflecting 
mirror 



Fully reflecting 
mirror 



Input Fields 



Output Fields 




Fig. 1, Schematic diagram of an OPA system. Here, the red solid lines 
represent the fields at the fundamental frequency and the blue dashed lines 
represent the fields at the second harmonic frequency. 

model for an OPA is as follows: 

8 = I, H x = 0, H 2 = i X (a* 2 a\ - {a\) 2 a 2 ) , 



L 



(40) 



e.g. see [20], [21]. Here ax and a 2 are the creation op- 
erators for the fundamental and second harmonic modes 
respectively. Also, kx > and k 2 > are parameters 
defining the reflectivity of the partially reflecting mirror at the 
fundamental and second harmonic frequencies respectively. 
Furthermore, % > is a parameter defining the strength of 
the x^ nonlinearity. 

The Hamiltonian in this model is a non-quadratic Hamil- 
tonian and so corresponds to nonlinear quantum system. The 
first term in the expression for H 2 can be interpreted as the 
annihilation of two photons at the fundamental frequency 
and the creation of a single photon at the second harmonic 
frequency. Similarly, the second term in the expression for 
H 2 can be interpreted as the annihilation of a single photon 



at the second harmonic frequency and the creation of two 
photons at the fundamental frequency. 

We will apply the theory developed in this paper to analyze 
the robust stability of this nonlinear quantum system. We 
first attempt to apply the results of Theorem [T] directly to 
this quantum system. Hence, we let 



M = 



and 



f(z, z*) = ix [zi (z\f - z\ (z|)J 



(41) 



where z\ — a\ and z 2 — a 2 . This defines a nonlinear 
quantum system of the form considered in Theorem [1] with 
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I. We now investigate whether this function /(•) 
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satisfies the conditions (0 and dTOb . First we calculate 



df(z,z*) 

dzx 
df(z,z*) 

dz 2 



d 2 f{z,z*) 
dz\ 
d 2 f(z,z*) _ 
dzi 



From this, we can immediately see that the conditions (0 
and ( fTOb will not be globally satisfied. 

In order to overcome this difficulty, we first note that 
any physical realization of a x^ optical nonlinearity will 
not be exactly described by the model (l40l but rather will 
exhibit some saturation of the nonlinear effect. In order to 
represent this effect, we could assume that the true function 
/(•) describing the Hamiltonian of the OPA is such that the 
first two non-zero terms in its Taylor series expansion ((3} 
correspond to the standard x^ Hamiltonian defined by fiTI ). 
Furthermore, we could assume that the true function /(■) is 
such that the conditions ^ and (TTOb are satisfied for suitable 
values of the constants 7 > 0, <5i > 0, 5 2 > 0. Here the 
quantity - will be proportional to the saturation limit. 

An alternative approach to dealing with the issue that 
the conditions (0 and ( TTOb will not be globally satisfied 
by the function /(•) defined in (HTt is to assume that 
these conditions only hold over some "domain of attraction" 
and then only conclude robust asymptotic stability within 
this domain of attraction. This approach requires a semi- 
classical interpretation of the function /(•) since formally 
the operators ax and a 2 are unbounded operators. However, it 
leads to results which are consistent with the known physical 
behavior of an OPA in that it can become unstable and 
oscillate if the magnitudes of the driving fields are too large; 
e.g., see [21]. In practice, the true physical situation will 
combine aspects of both solutions which we have mentioned 
but we will concentrate on the second approach involving a 
semi-classical "domain of attraction". 

In order to calculate the region on which our theory can 
be applied, we note that for our OPA model, the condition 



(|9]i will be satisfied if 
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We now consider the case in which ||zi| 
case, the condition d42l is equivalent to the condition 
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In the case that II 



4 7 2 X 2 



the left hand side of d42l) 
is always negative and the right hand side of (l42l is always 
positive. Hence in this case, the condition d42b will always 
be satisfied. Also, the condition ( fTOb will be satisfied if 



(44) 



The conditions d43l and d44t define the region to which 
our theory can be applied in guaranteeing the robust mean 
square stability of the OPA system. This region is represented 
diagrammatically in Figure [2] The constraints (|43> and (l44l 
can be interpreted as bounds on the average values of the 
internal cavity fields for which robust mean square stability 
can be guaranteed; see also [22]-[24]. 




Fig. 2. Diagrammatic representation of the region for which the constraints 
(9) and )10t are satisfied. Here A = „ + a / Tzkzz + ^1 ■ 
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We now investigate the strict bounded real conditions (fT3l . 
For this system, it follows from the definition (fT3l that 
the matrix F is given by 
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which is Hurwitz for all k± > and «2 > 0. Thus, the 
condition ( fT3l is always satisfied. Also, we calculate the 
transfer function matrix E (si — F) 1 as 



Eisl-Fy 1 J& 



6+^ 







-1 



It is straightforward to show that this transfer function matrix 
has an H°° norm of 



E(sl - FY 1 JE ] 



Thus, for this system, the condition 
condition 





/ max < 


r 4 












00 




L«? 





is equivalent to the 



< 



(45) 



Hence, using Theorem Q] and Lemma [5] we can conclude 
that the OPA system d40b is robustly means square stable 
provided that the condition (|43T > is satisfied and Heisenberg 
evolution of the quantities ai(t) and 02 (t) are such that the 
conditions d43l and d44l i remain satisfied. 

Note that in most experimental situations, K2 > «i; e.g., 
see [22]. This means that 



4 4 



2 



If we then equate ^_ — 2 5 we obtain 7 = 7^- which can be 
substituted into the right hand side of d43l to obtain an upper 
bound on the region for which the conditions ((9) and (TToT > are 
satisfied. This region is defined by d44l i and the inequality 



-2 



< 



Si 



IN1II 2 



16 X 2 



4|ki!l 2 -ii2 



(46) 



Also, note that the region defined by (l43l and (|46] > will only 
be an upper bound on a domain of attraction for the OPA 
system. To find an actual domain of attraction for this system, 
we would need to find an invariant subset contained in the 
region defined by (1431 and d46l l. Such an invariant set could 
be chosen to be an ellipsoidal region defined by the quadratic 
Lyapunov function arising from the matrix P solving d29| i. 

IV. Conclusions 

In this paper, we have extended the robust stability result 
of [15] to the case of non-quadratic perturbations to the 
Hamiltonian which depend on multiple parameters. This led 
to a robust stability condition of the form of a multi-variable 
small gain condition. This condition was then applied the 
robust stability analysis of a nonlinear quantum system 
consisting of an OPA and the stability region for this system 
was investigated. 
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